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ABSTRACT 



The deformation of an unstable baroclinic wave as a result of mean 
horizontal wind shear is studied with the use of a five- layer primitive 
equation model. Two basic wind profiles were tested: 1) uniform zonal 

flow, 2) jet type zonal flow. The jet profile produced a cold front 
and a warm front, with the warm front being the more intense. The 
uniform zonal flow produced a cold front. The cold fronts were pro- 
duced by horizontal shearing deformation fields and the warm front was 
produced by horizontal stretching deformation fields. The formation of 
the warm front was also enhanced by strong horizontal convergence. The 
occlusion process was not observed in the experiments. 
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I , INTRODUCTION 



The concept of atmospheric fronts has been accepted for many years, 
yet at this time there is no mathematical theory which can explain the 
formation of a front. Williams [1967] has indicated that the basic 
problem is finding a process which can produce a small-scale motion 
field from one which initially contains only large-scale motions. 

A number of studies have been conducted evaluating mechanisms 

which are contributing factors in the formation of a front from an 

initial large-scale baroclinic wave. Stretching deformation wind 
* 

fields have been studied by Stone [1966] and Williams and Plotkin 
[1968], while shearing deformation wind fields have been studied by 
Arakawa [1962] and Williams [1967]. Stone [1969] indicates that these 
effects should be of about the same importance. Arakawa [1962], 
Edelmann [1963], Williams [1965, 1967] and 0kland [1969] have studied 
atmospheric frontogenes is with vertical deformation fields with the 
general result that a realistic front is produced within a reasonable 
time period. Bushby and Timpson [1967] and Bushby [1969] have studied 
the dynamics of fronts and frontal rainfall with a primitive equation 
model utilizing real data with good results in a twenty-four hour 
forecast of rainfall. 

The purpose of this study is to investigate the deformation of a 
unstable baroclinic wave as a result of horizontal shear in the mean 
wind field. The primitive equation model is basically the model used 
by Harrison [1969] with a slightly modified lower boundary condition 
and different geographical and theoretical areas of application. 
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II. THE PRIMITIVE EQUATIONS 



Consider the primitive equations in the following form, with 
heating, friction and curvature terms neglected and with the Coriolis 
parameter a linear function of y: 



Su 
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where 

L( S ) = ^ (us) + (vs) + (COS), 
and s = u, v, 9. 




Equations (1) through (3) and equation (6) are prognostic 
equations. Equations (1) and (2) are the momentum equations, 
equation (3) is the first law of thermodynamics and equation (6) 
is the lower boundary condition. The diagnostic equations (4) and 
(5) are the mass continuity and hydrostatic equations, respectively. 

Equation (6), the lower boundary condition was obtained by 
setting d0/dt = gw = 0, which implies that w = 0 at the lowest level. 



12 



Phillips [1963] has shown that the error resulting from the neglect of 
vertical motion at the lowest level is very small. 
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III. FINITE DIFFERENCES 



Equations (1) through (6) were solved numerically by introducing 
finite differences in x, y, p and t. The space finite differences 
employed were a variation of a scheme constructed by Arakawa [1966] 
[also see Lilly, 1965]. The basic scheme is designed to conserve 
total energy and also to conserve the squares of u, v and 9 with 
respect to the advection terms. 

The vertical domain is divided into five layers with the arrange 
ment of variables as shown in Fig. 1. 
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Fig. 1. Distribution of variables in the vertical. 

The finite difference form of the L(s) operator can be written 



L(s) i,j,k 



L^ U i+l,j,k + U i,j,k^ ^ S i+l,j,k + S i,j,k^ 

^ U i,j,k + U i-l,j,k^ ^ S i,j,k + S i-1, j .k^W^* 
+ L^ V i,j+l,k + V i,j,k^ ^ S i,j+l,k + S i,j,k^ 

■ <v t,j,k + v i, J -i,k ) < s i, J ,k + s i,j-i,k ) ] /4i >’ 
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The linear space derivatives are of the general form 
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The hydrostatic equation (6) can be written in finite difference 
form following Lorenz [1960] as 

(9. . . + 9. . . ..) 

0 = 0 + LtrJs+g- 

i,j,k+2 1 °i,j,k^ 2 



C P [£)" • &$}■ 



The form of vertical finite differences as specified in the L(s) 

i * * 

operator and the hydrostatic equation were modified at level zero in 
the following manner. 

1. 9, u and v were linearly extrapolated from levels one and 
three to level zero. 

2. 50/ 5p at level zero was approximated by a one-sided 



difference . 

Centered time differences were used throughout this study with the 
exception of the first time step which is forward. To prevent sepa- 
ration of solutions, the program was restarted every twelve hours with 
a forward time step. 
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The following basic steps were used in integrating the given 
system of equations. 

1. The potential temperature and wind fields at the odd- 
numbered levels and the geopotential of the 1000 mb surface are 
prescribed initially . 

2. The mass continuity equation is integrated in space from 
the highest level to the surface to obtain od. 

3. The hydrostatic equation is integrated in space from the 
surface to the highest level to obtain 0. 

4. The momentum equations are integrated in time to obtain 
u and v at time t + At . 

5. The first law of thermodynamics is integrated in time to 
obtain 0 at time t + At . 

6. The equation for the lower boundary condition is integrated 
in time to obtain the 1000 mb geopotential at time t + At . 

7. Steps 2 through 6 are repeated for the number of time steps 

desired . 
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IV. BOUNDARY CONDITIONS 



The boundary conditions are specified in a manner that conserves 
mass and energy within the grid domain. At the northern and southern 
boundaries the normal component of velocity is zero and periodicity is 
enforced at the eastern and western boundaries. Fig. 2 illustrates the 
horizontal grid domain. 
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Fig. 2. Horizontal grid 

The east-west boundary condition forces cyclic continuity on all 
parameters. This is accomplished by requiring the values in the two 
columns on the right to be the same as the values in the two columns 
on the left. This condition is specified by 

1>J M-l,j 2 , j M , j 

where 



s = u,v,0 5 0. 

The north-south boundary condition allows flow along the boundary 
but does not allow flux across the boundary. This is accomplished by 
imposing the following condition on the v term of the advection operator, 



or 



v + v = 0 

i,l i,2 



V i,l ~ ’ V i>2 



V. + v. - 0 

i,N i,N-l 



v — - v 
i,N i,N-l 
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The values of geopotential and potential temperature in the outer row 
are set equal to the value in the next interior row. This is speci- 
fied by 



where 




S i,N S i,N-l 



(7) 



s = 0,9. 

As a result of v being set equal to zero between the two outer grid 
points along the northern and southern boundaries, values of u and O) 
are not required outside the boundaries. 

The upper boundary condition of cd = 0 eliminates flux through the 
top of the domain. The lower boundary condition has been specified in 
equation (6). 

Various boundary conditions were tested (see Appendix A) before 
acceptance of the above method. 
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V . INITIAL CONDITIONS 



The initial fields of u, v and 0 were specified in order to 
generate a baroclinic zonal current. A small amplitude disturbance 
in the geopotential field; which is independent of pressure, was 
superimposed on the basic zonal flow. The mean zonal current (u) has 
the following forms; 

1 ) jet profile 



u = 1; (p ■ p 0 ) sin2 OF) 



2 ) uniform zonal profile 

- du , . 

u = — (p - p ) 



bp 



A 



where bu/bp is constant. The initial data fields were generated in 
the following steps. 

1. The geopotential at p = p was calculated from 



0 = f U ~~ cos (*&) sin^ ( 77 ^) 



2tt 



\ L / 



VW J 



( 8 ) 



2. The potential temperature fields were specified from 



0 = 



P 

o\ 






R 



|j + 9. + | (p - pj. 



where d 0 /c)p is a constant, 



For uniform zonal flow 



J = f o y + 3 o y /2. 
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For jet flow 



J = f Q y + 5y 2 /2 




3. The hydrostatic equation was then solved for the geopotential 
at all required levels. 

4. The initial wind fields were then calculated from the geo- 
potential fields using the geostrophic approximation. Thus, the 
initial fields were in geostrophic balance at time zero. 

The basic current is zero at the lowest level and increases at a 
rate of 5 mps/100 mb to a maximum of 45 mps at 100 mb. Therefore 
cki/dz decreases with increasing altitude which reduces the baroclinic 
instability in the upper layers. This approximates the effects of a 
tropopause on the dynamics. 

The initial disturbance amplitude was 10 meters or .125 mb. The 
small initial amplitude was chosen to prevent the excitation of large 
amplitude gravity waves that occur since the initial wind field was 
geostrophic . 

The value of d0/dp used was -.376°C/mb. This value is approxi- 
mately one-half the value determined by Gates [1960] for winter over 
the United States, but it has the advantage of forcing the fronts to 
develop more rapidly and thus to reduce the computer time required. 

The value used approximated the standard atmosphere quite well in the 
lower levels, see Table I, but is less stable than the standard atmos- 
phere in the upper levels. The stratosphere is clearly not represented 
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by this temperature field, however the dynamical presence of the 
stratosphere is represented by u. 



Pressure (mb) 


Height (meters) 


Jet 


Standard 

Atmosphere 


100 


15325 


16210 


300 


9112 


9160 


500 


5576 


5570 


700 


3015 


3010 


900 


988 


990 


1000 


110 


110 



Table I. Average initial height of constant pressure surface. 

The beta plane approximation was made, where the coriolis parameter 
is a linear function of y. The quantities f^ and (3^ were computed at 
45° north latitude. 

Various combinations of grid size, time step and wind profile were 
tested. Table II is a summation of these various experiments. 
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Table II. Summary of experiments. 



VI. RESULTS 



All experiments were conducted using the values of du/dp, S0/^p 
and initial disturbance amplitude as specified by equation (8) in the 
initial condition section of this paper. Table II gives a listing of 
the specific variations in each of the following experiments. 

Experiment No. 1 . This is the basic jet flow experiment with a 
grid size of 200 km in x and y. The initial 1000 mb height analysis 
is given in Fig. 3 and the initial height and potential temperature 
analysis for the 900 mb and 700 mb surfaces are given in Figs. 4 and 
5. The growth of the systems from day one through day four at the 900 
mb surface is given in Figs. 6 through 9. By the end of three days a 
cold front and a warm front have formed and the temperature gradients 
are large by the fourth day. Note that the warm front contains larger 
gradients than the cold front. Large convergence (not shown) is associ- 
ated with the warm front. The 700 mb analysis of height and potential 
temperature for days three and four is given in Figs. 10 and 11. The 
cold and warm fronts have the proper tilt in the vertical, with a 
slope between 1/125 and 1/250, which is computed at day three and for 
both fronts. The growth of the surface pressure centers, divergence 
and vorticity at 900 mb is given in Table III. It should be noted in 
Table III the asymmetry between the maximum and minimum values of each 
quantity . 

Experiment No. 2 . This is the basic uniform zonal flow experiment 
with a grid size of 200 km in x and y. The initial 1000 mb height 
analysis is given by Fig. 3 and the initial height and potential 
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mb 

rs (m) 


Divergence 
(x IQ-^ti 


Vorticity 
(x 10' 5 ’) 


liAI 


Low 


High 


Max 


J Min 


Max 


Min 


Initial : 


100 


120 


, 3 


-.3 


.7 


-.7 


i 


90 


130 


2,2 


-2.4 


1.0 


-1.0 


2 


40 


160 


3.7 


-7.3 


2.9 


-1.8 


3 


-60 


190 


6.5 


-17.1 


9.4 


-3.0 


4 


-200 


270 


20.7 


-41.9 


18.9 


-6. 6 



Table III. Experiment No. 1* 

temperature fields for the 900 mb and 700 mb surfaces are given in 
Figs. 12 and 13. The growth of the systems from day one through day 
four at the 900 mb surface is given in Figs. 14 through 17. By day 
three a cold front has formed to the west of an apparent warm front. 
However there is a single vorticity maximum in the region and there 
is no evidence of a warm front at higher levels. The front appears 
to be an occlusion but it did not evolve from two separate fronts with 
associated vorticity maxima. It was also noted that in the northern 
region of the low there was a packing of potential temperature lines 
which somewhat resembles the warm front in the jet flow experiment. 

The 700 mb pressure analysis of height and potential temperature for 
days three and four are given in Figs. 18 and 19. The front tilts to 
the west in the vertical with a slope between 1/60 and 1/120 at day 
three. The growth of the surface pressure centers, divergence and 
vorticity at 900 mb is given in Table IV„ 

Experiment No. 3 . This experiment is the same as experiment 1 
except that the grid length was reduced to 100 km in x and y. The 
results are basically the same as in experiment 1, with the exception 
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8.3 


-8.8 
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10.0 


-2.9 


4 
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340 


23.3 


-43.7 
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-4.5 



Table IV. Experiment No. 2. 

that the growth of the systems is now slightly faster with a slightly 
larger vorticity. The 900 mb height and potential temperature analysis 
for days three and four are given in Figs. 20 and 21. Note that the 
temperature gradients are larger in this experiment at four days than 
in experiment 1. This indicates that the larger truncation error with 
the larger grid size inhibits the energy cascade to small scales. 

Experiment No, 4 . This experiment is the same as experiment 2 

* 

except that the grid size was reduced to 100 km in x and y. A front 
similar to the one in experiment 2 developed, but in addition, a 
gravity wave with a wavelength of approximately 188 km was super- 
imposed on the height analysis. The gravity wave was initially 
observed at day two close to the north and south boundaries, and by 
day three had completely covered the domain. It was apparently excited 
by the large shear at the boundary as generated by equation (7). 

Experiment No. 5 . This experiment is the same as experiment 1 
except that Ay was reduced to 100 km and W was reduced to 2000 km. 

Cold and warm fronts developed within three days that were similar to 
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the ones generated in experiment 1. Superimposed on the long wave 
pattern was a very small amplitude gravity wave with a wavelength of 
approximately 188 km. 

Experiment No. 6 . This experiment is the same as experiment 2 
except that Ay was reduced to 100 km and W was reduced to 2000 km. 

The results are the same as in experiment 2, but in addition, a 
gravity wave of small amplitude and approximate wavelength of 188 km 
was superimposed on the long wave pattern. 

Experiment No. 7 . This experiment is the same as experiment 1 
except that Ay was increased to 400 km and W was increased to 8000 km. 
See Figs. 22 and 23 for day three 900 mb height and potential tempera- 
ture analysis. By day three a cold and warm front, similar to the 
results of experiment 1, had formed with the disturbance being limited 
to the central region of the domain. 

Experiment No. 8 . This experiment is the same as experiment 2 
except that Ay was increased to 400 km and W was increased to 8000 km. 
Gravity wave propagation, with an approximate wavelength of 400 km, 
from the north and south boundaries was too strong to permit an 
evaluation of the long wave system. 

In none of the experiments was there an indication of a front at 
500 mb or above. All energy calculations (see Appendix B) indicated 
that the potential energy was three orders of magnitude greater than 
the kinetic energy and that there was a small energy sink in the 
system up to four days. 
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Minimum contour 100 meters, contour interval 
10 meters. 



Fig. 3. Initial 1000 mb height analysis. 
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Fig. 6. Experiment No. 1, 24 Hour 900 mb height and potential temperature analysis. 
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Fig. 7. Experiment No. 1, 48 Hour 900 mb height and potential temperature analysis. 
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Fig. 8. Experiment No. 1, 74 Hour 900 mb height and potential temperature analysis. 
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Fig. 9. Experiment No. 1, 96 Hour 900 mb height and potential temperature analysis. 
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Fig. 11. Experiment No. 1, 96 Hour 700 mb height and potential temperature analysis. 
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Fig. 13. Experiment No. 2, Initial 700 mb height and potential temperature analysis. 
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Fig. 14. Experiment No. 2, 24 Hour 900 mb height and potential temperature analysis. 
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Fig. 16. Experiment No. 2, 72 Hour 900 mb height and potential temperature analysis. 
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Fig. 18. Experiment No. 2, 72 Hour 700 mb height and potential temperature analysis. 
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Fig. 20. Experiment No. 3, 72 Hour 900 mb height and potential temperature analysis. 
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Minimum contour 800 meters, interval 30 meters. 



Fig. 22. Experiment No. 7, 72 Hour 900 mb height analysis. 
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Fig. 23. Experiment No. 7, 72 Hour 900 mb potential temperature 
analysis . 
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VII. CONCLUSIONS 



The time evolution of a baroclinicly unstable disturbance was 
studied for two basic zonal currents. With horizontal wind shear, 
cold and warm fronts developed and are similar to the results of 
Edelmann [1963] and 0kland [1969]. With no horizontal wind shear, 
a cold front formed which was similar to the results of Williams 
[1969]. This cold front resembled an occluded front at the surface 
but did not have the life history of an occlusion. In general, the 
warm -front formed from the stretching deformation discussed by Stone 
[1966] while the cold fronts formed from the shearing deformation 
discussed by Williams [1967]. In the jet case, the warm front was 
further enhanced by strong convergence through the action of the 
vertical deformation field. The lack of frontogenes is at upper levels 
is explained by the reduction of the vertical wind shear at those 
levels. The classical occlusion process was not observed in any of 
the experiments. 

Further studies with this model should be conducted with the 
addition of vertical and horizontal diffusion and with the addition 
of latent heat of condensation. Other basic currents should be 
tested, including those with an explicit tropopause . 
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APPENDIX A 



EXPERIMENTS WITH BOUNDARY CONDITIONS 

Before acceptance of the north-south boundary conditions as 
specified in Chapter V, various other schemes were tested with 
limited success. A sample of these schemes with the results follows. 

1. The boundary was placed at the outer grid points with the 
values of all variables determined from a linear extrapolation from 
the next two interior grid points. The integration became unstable 
after approximately 120 six -minute time steps. 

2. The boundary was again placed at the outer grid points and a 
special one-sided difference modification was made to the basic diag- 
nostic and prognostic equations. Strong unrealistic wind shear 
developed along the boundary and the numerical integrations eventually 
became unstable. 

3. Cyclic boundary conditions were applied in the north-south 
as well as in the east-west directions. The discontinuity in the 
geopotential field forced the use of a complete wave in both the 
north-south and east-west directions. This doubled the number of grid 
points in the domain and made the size of the computer program too 
large for practical use. 
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APPENDIX B 



ENERGY COMPUTATIONS 

In pressure coordinates, the total potential energy may be 
written as 



EP 



■i If 



C T dp dA, 
P 



B.l 



A o 



where n dA M represents an area increment. Now, substituting the 
definition of potential temperature in (B.l), the expression for 
the total potential energy becomes, 



EP 



■ g 1 i C P e Or) dA - 



B.2 



A o 



In pressure coordinates, the total kinetic energy is written 



as 



EK 



= h f f (u2 + y2) dp dA * 



B.3 



A o 
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